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Ground state properties of the Kondo model for manganese oxides in one dimension are studied 
using numerical techniques. The large Hund coupling (Jh) limit is specially analyzed. A robust 
region of fully saturated ferromagnetism (FM) is identified at all densities. For open boundary 
conditions it is shown exactly that the ground state is FM at Jh = oo. Hole-spin phase separation 
competing with FM was also observed when a large exchange J between the Mn 3+ ions is used. 
As the spin of the transition metal ion grows, the hole mobility decreases providing a tentative 
explanation for the differences between Cu-oxides and Mn-oxides. 

PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 
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Doping of transition metal oxides with perovskite 
structure induces remarkable phenomena such as high 
temperature superconductivity, charge ordering, and 
anomalous transport properties. Typical examples are 
the layered Cu-oxides and Ni-oxides. Recently, an- 
other dramatic property of doped carriers in transition 
metal oxides has been revealed: at low temperature 
(T) La\- x C a x MnOz changes from an antiferromagnetic 
(AF) insulator to a ferromagnetic metal as x grows from 
0. "Colossal" magnetoresistance phenomena is observed 
as T increases from the ferromagnetic phase At 
the critical ferromagnetic temperature, Tf M , a metal- 
insulator transition occurs and the high temperature 
regime is a paramagnetic insulator. At larger doping 
x > 0.5, a charge-ordered AF state was detected 0. 

The double-exchange (DE) mechanism has been used 
to explain the FM phase in Mn-oxides ||, since the mo- 
bility of the e g electrons increases for a polarized back- 
ground of localized ti g electrons. However, recent reex- 
amination of the DE model, where electrons move with 
an effective hopping proportional to the cosine of half the 
angle between the localized spins (assumed classical), re- 
vealed that electrons can actually collect a Berry phase 
that may induce states of low energy in the spectrum [Q . 
This effect could reduce the Tf M and alter the behav- 
ior of the system at higher temperatures. In addition, 
from recent numerical studies of the Kondo Hamiltonian 
for Mn-oxides it was concluded that the ground state at 
nonzero hole density is not FM but a spin singlet, al- 
though with local ferromagnetic correlations ||. These 
results suggest that electronic models for the Mn-oxides 
more realistic than the DE model may have a phase di- 
agram richer than expected which deserves to be stud- 
ied with state-of-the-art computational techniques. Such 
studies would clarify, among other issues, whether FM 
indeed exists in the T = ground state and what phases 
are in competition with it. In addition, it would be im- 
portant to analyze the origin of the drastic differences 
between Cu-oxides and Mn-oxides regarding their doped- 



induced properties. The cuprates have metallic and su- 
perconducting phases, while the manganites are either 
charge-ordered insulators or ferromagnets [Q . In this pa- 
per, both issues are addressed. 

The Hamiltonian widely used for Mn-oxides is p] 
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Here the first term represents the e g -electron transfer be- 
tween nearest-neighbor Mn-ions at sites m and n. The 
second term corresponds to the Hund coupling ( Jh > 0) 
between the S = 3/2 t2 g localized spin S n and the spin 
<7 n of the mobile e s -electron at the same site. Coulom- 
bic repulsion in the e g -band is not included in Eq.(l) 
but it is considered in some calculations below ||. Phc- 
nomenologically, it is expected that Jh 3> t which favors 
the alignment of the itinerant and localized spins. For 
Mn 3+ , where the e g level is occupied, the resulting spin 
is 2, while for Mn 4+ (vacant e g state) the spin is 3/2. 
Thus, at Jh = oo the effective degrees of freedom become 
S — 2 "spins" and ,5 = 3/2 "holes" which is the language 
convention followed below. In this limit Eq.(l) reduces to 
a simpler hole-hopping Hamiltonian. For "spin" values 
S > 1, the hopping term at Jh — oo is M 



Hj H =oo = -ty] P mn Qs(y), 



(2) 



(mn) 



where y = S m • S n /5(5-l/2) and for Mn-oxides (5 = 2) 
the polynomial Qs[y) is 
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For Ni-oxides (S = 1) Q x (y) = (1 +y)/2 §. Hamil- 
tonian Eq.(2) is rotational invariant and it acts only on 
links containing one spin and one hole. P mn permutes 
the hole and the spin. The S n operators are standard 
and in Eqs.(2,3) they act over both the spin and hole in- 
volved in the hopping process, i.e. the kinetic term not 
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only interchange their sitepositions but it can flip their 
spin projections as well y,|7|. This is an important differ- 
ence with respect to the standard t — J model. At finite 
but large Jh, Eq.(2) is supplemented by a Heisenberg 
interaction between nearest neighbor spins of the form 
J^{ mn )(Sm ' S n — n m n n /A), in the standard notation. 
H This term only affects the S spins since in hole states 
the orbital with mobile electrons (e g ) is empty. 

As numerical techniques the exact diagonalization ||] 
(ED) and the infinite and finite-size density matrix renor- 
malization group || (DMRG) methods were used. The 
special case 5 = 1 will be studied together with the rel- 
evant case for Mn-oxides (5 = 2) to observe the effect of 
different transition metal spins in the phase diagrams. In 
addition, experience with Cu-oxides models shows that 
the analysis of Hamiltonians in ID chains provides insight 
on qualitative features of the phase diagram that survive 
the increase in dimensionality. Since ID is the best for 
numerical techniques, our study below is restricted to 
chains |u|. Note that the special case 5=1 and ID has 
intrinsic relevance since the compound Yi^ x Ca x BaN iO^ 
with NiO chains has been recently synthesized jy] . 

First let us show that the ground state of Eqs.(l,2) in 
ID is a fully saturated ferromagnet (FM) for hole density 
x 7^ 0, 1 in the case of open boundary conditions (OBC) 
and Jh — oo (or J = 0) . The proof is a simple extension 
of Kubo's results for the DE model E3j: the matrix ele- 
ments of the J = Hamiltonian Eq.(2) in the standard 
5 2 -basis are nonpositive, which is shown by explicitly 
writing the polynomials Q s (y) in matrix form for a given 
link. For a subspace with a fixed total spin projection 
S total > th^ s property implies that the ground state is non- 
degenerate and the coefficients of this state expanded in 
the 5 2: -basis are of the same sign(Perron-Frobenius the- 
orem) Q. Since in each subspace with a fixed S^ otal 
the state \S max , $ total) with the maximum possible spin 
S m ax has all of the coefficients of the same sign, then the 
ground state in such subspace must be the ferromagnetic 
state. Q.E.D. As in Kubo's results, this proof is not valid 
for twisted boundary conditions, including periodic and 
antiperiodic conditions (PBC and APBC, respectively), 
since in this case across the boundary it is possible to 
move fermions from site N to 1 (for a N-site chain) col- 
lecting anticommutation signs which make the matrix el- 
ements not necessarily negative. For the same reason 
the proof is not extended to higher dimensions even with 
OBC. Another important detail is that each S^ otal sub- 
space is assumed to be fully connected (i.e. all elements 
of the basis are reached by successive applications of the 
Hamiltonian over any of its members). Actually in the 
ID t — J model, the Hilbert space is disconnected for 
J = and the theorem does not hold. However, if 5 > 1 
the Hilbert space is connected since together with a hop- 
ping, the spins and holes can flip their spin projections 

Solving exactly Eq.(2) on finite chains with OBC we 
verified the existence of ground state FM for any x ^ 



0, 1. We also observed that using PBC(APBC) for an 
odd(even) number of holes (closed shell) for both 5=1 
and 2, the ground state at J = is again a FM. Thus, to 
increase the size of the chains that can be studied numer- 
ically with ED techniques, preserving the correct proper- 
ties of the ground state, below closed-shell BC have been 
used in addition to OBC. On the other hand, using PBC 
for an even number of holes, the ground state is a spin 
singlet Thus, in principle these BC should not be 
used for finite cluster studies of T = properties p6[ . 
However, the boundary dependence of the ground state 
total spin on finite chains suggests the existence of low 
energy singlet excitations close to the ground state [|],|| . 

In Fig. la the phase diagram for 5 = 1 is shown. Chains 
with up to 16 sites were studied with ED. Using the 
closed-shell BC, a robust region of FM was identified. 
The FM boundary was also calculated using the DMRG 
technique at x = 0.25, 0.50 and 0.75 with chains of up to 
40 sites and OBC . The critical coupling to reach FM 
obtained with ED and DMRG are in excellent agreement, 
showing that size effects are small |1| . 

Also in Fig. la, regions of phase separation (PS) and 
hole binding (B) have been identified using ED by cal- 
culating the compressibility and binding energy from 
E(M ± 2), E{M ± 1), and £(M),where E(M) is the 
ground state energy for M holes |Q. PS is character- 
istic of spin- hole models, although not contained in the 
DE model, and it occurs at large J/t when the kinetic 
energy gained by making holes mobile is overcome by the 
magnetic energy lost. This has been observed in the t — J 
model in both ID and 2D || and it is not surprising to 
find a similar behavior in models with a larger 5. How- 
ever, at x = 0.5 PS appears at J ~ 1.1 which is much 
smaller than the coupling J/t ~ 3.0 needed in the t — J 
model p0[ , suggesting that electronic PS may play an 
important role in the physics of Mn-oxides. Close to PS 
there is a region (B) where holes form mobile pairs. Su- 
perconductivity is likely in this regime which also appears 
in the t — J model [^lj . 

In Fig. lb, the 5 = 2 phase diagram is shown. Here the 
chains accessible to ED have 8 sites, while with DMRG 
we studied up to 30 sites. The finite-size effects on the 
FM line are small. As for 5 = 1, a robust regime of 
FM has been identified. We observed that the FM lines 
of phase transitions of Fig. la, b are very similar if J is 
scaled as J5(5 + 1). PS and B are also present. Note 
that for 5 = 2 the intermediate metallic regime is very 
narrow, i.e. as 5 grows the trend observed in Fig. la, b 
suggests that the PS and FM phases become dominant. 
To understand this effect note that the "vaporization" of 
a large spin cluster in the PS region, held together by 
magnetic forces of effective strength J5(5 4- 1), can only 
occur if a substantial gain in kinetic energy, dominated 
by the scale t, is achieved by the vaporized spins in the 
hole rich region. Thus, the critical J/t\ps leading to PS 
should follow J/t\ps ~ 1/5(5 + 1), or a faster decay with 
5 once the reduction in hole kinetic energy due to the 
scattering with the magnetic background is considered. 
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This argument breaks down at very small J/t where the 
FM state is energetically competitive. Thus, in the large 
S limit the metallic regime at T = would disappear in 
favor of PS and FM phases. 

The study of the intermediate metallic regime is in- 
teresting since it may contain features of the param- 
agnetic insulator observed in manganite experiments at 
T > T C FM . As S grows, a tendency towards localiza- 
tion is indeed observed in this regime. For example, in 
Fig. 2a the bandwidth W of one hole for S = 1/2, 1, and 
2 vs t/JS(S + 1) is shown. It is clear that W rapidly 
diminishes with S. Similar effects are observed in Fig. 2b 
where the kinetic energy in the metallic region is pre- 
sented at x = 0.5. Such a mobility reduction is rea- 
sonable since holes destroy the magnetic order in their 
movement, paying an energy that rapidly grows with 
the spin. The large effective mass for S = 2 indicates 
that small perturbations away from a translational in- 
variant system may localize the carriers. This effect and 
the dominance of PS and FM at large S is a possible 
explanation of the notorious experimental differences at 
low-T between Cu-oxides, which show metallic and su- 
perconducting phases, and Mn-oxides, which have FM 
and charge-ordered phases H (note that previous studies 
in Ni-oxides have shown that char ge-o rdering and phase 
separation are related phenomena |22| ). 

The Hamiltonian Eq.(l) does not include a spin-spin 
exchange between the ti g spins since such a coupling 
arises indirectly from Eq.(l) mediated by the e g elec- 
trons. However, as x — > 1 this RKKY coupling vanishes, 
which is not in agreement with experiments that show AF 
order at x = 1. To remedy this problem phenomenologi- 
cally, a Heisenberg interaction with coupling J' between 
the "holes" (which carry spin) can be introduced. At 
x = 1 now the ground state would be AF as for x = 0. 
For completeness, the numerical analysis of Fig. la was 
repeated adding a term J' J2(mn) (Sm • S n — n m n n /4) be- 
tween the 5 = 1/2 holes and considering J' = J (see 
Fig. 3a). Up to x » 0.5 the phase diagram is similar to 
the results of Fig. la, but as x grows further the coupling 
needed to reach the FM phase is reduced, as expected. 
FM remains robust at all densities and small J/t. 

Finally, in Fig. 3b the phase diagram of Hamiltonian 
Eq.(l) solved exactly on a 4-site chain and using DMRG 
on an 8-site chain, both with OBC and at x = 0.5, is 
shown. Here, an on-site Hubbard repulsion U ft for e g 
electrons is incorporated in the Hamiltonian. In agree- 
ment with our analysis of the Jh — oo limit, the ground 
state is FM at large Jh- As U/t grows, FM becomes 
stable even at small Jh which is natural since in the 
Hubbard model ferromagnetism tends to be favored in 
strong coupling. As with model Eq.(2), the non-FM re- 
gion contains a small ferromagnetic component, i.e. for 
Jh = U = 1 and N — 8 using DMRG we find spin one 
in the ground state. 

Summarizing, a numerical study of the ID Kondo 
Hamiltonian for Mn-oxides in the Jh ^ t limit and also 



at finite Jh has been presented. The phase diagram con- 
tains a robust FM region at all densities in agreement 
with expectations from the DE model, and contrary to 
the case of the ID t — J model for Cu-oxides. For OBC 
we showed exactly that the ground state at Jh = oo and 
x 0,1 is FM with maximum spin. At finite exchange 
J, a narrow metallic window with low mobility carriers 
separates the FM region from a phase separated regime, 
not contained in the DE model. As the spin S grows, 
PS and FM become dominant, the intermediate metallic 
regime tends to disappear, and the mobility of holes in 
the non-FM region rapidly decreases. 
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Figure Captions 

1. (a) Phase diagram of the 5 = 1 Hamiltonian Eq.(2) 
on a ID chain (JS(S + 1) is the coupling scaled 
with the spin, x the hole density, and t = 1). The 
gray line and open points denote the FM phase 
boundary. Open circles(lO), rotated triangles(lO), 
squares(12), diamonds(14), and triangles(16) cor- 
respond to ED results using closed-shell BC (in 
parenthesis the number of sites). DMRG results for 
the FM boundary with OBC and N = 40 are also 
shown (stars). The solid line joining full circles is 
the PS boundary calculated with ED (N = 12). A 
hole binding (B) region exists between the dashed 
line and PS; (b) Same as (a) but for the 5 = 2 
Hamiltonian. The boundaries of PS and FM were 
evaluated with ED on chains of N — 8 sites. 
DMRG results on N = 30 OBC are shown (stars). 
In the intermediate metallic region there is binding. 



2. (a) Bandwidth W/J vs t/JS(S + 1) for one hole 
and 5 = 1/2 (full diamonds), 1 (open circles), and 
2 (full squares) calculated using ED techniques on 
N = 12, 12, and 8 chains, respectively; (b) Ki- 
netic energy (hopping term ground state expecta- 
tion value) in units of t vs t/JS(S+ 1) obtained at 
x = 0.5 for 5 = 1/2, 1, and 2 using N = 12, 12, and 
8 chains, respectively. 

3. (a) Same as Fig. la but adding a Heisenberg in- 
teraction between the 5 = 1/2 holes of strength 
J' = J; (b) Phase diagram of the Kondo model 
Eq.(l) at x = 0.5 including an e g Hubbard repul- 
sion of strength U (t = 1). Open (full) squares are 
results obtained with ED (DMRG) on chains with 
N=4 (8) sites. The non-FM region indicates weak 
ferromagnetism in the ground state, according to 
the DMRG results. 
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